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Abstract 

We systematically analyze quantum corrections in see-saw scenarios, including effects 
from above as well as below the see-saw scales. We derive approximate renormaliza- 
tion group equations for neutrino masses, lepton mixings and CP phases, yielding an 
analytic understanding and a simple estimate of the size of the effects. Even for hier- 
archical masses, they often exceed the precision of future experiments. Furthermore, 
we provide a software package allowing for a convenient numerical renormalization 
group analysis, with heavy singlets being integrated out successively at their mass 
thresholds. We also discuss applications to model building and related topics. 
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1 Introduction 



The observed smallness of neutrino masses finds an attractive explanation in the see-saw 
mechanism PQ El IH IS] • The hght neutrino masses are, at tree-level, given by the famous 
see-saw relation 

= -{ni^'''^f M-^ m^^''^ . (1) 

This relation emerges from integrating out heavy, singlet neutrinos with mass matrix 
M. The Dirac neutrino mass m^^'^^'^ is proportional to the neutrino Yukawa coupling 
Yy. Clearly, the see-saw operates at high energy scales while its implications are measured 
by experiments at low scales. Therefore, the neutrino masses given by Eq. (P) are subject 
to quantum corrections, i.e. they are modified by renormalization group (RG) running. 

The running of neutrino masses and lepton mixing angles has been investigated inten- 
sively in the literature. For non-hierarchical neutrino mass spectra, RG effects can be very 
large and they can have interesting implications for model building. For example, lepton 
mixing angles can be magnified [HI 13 El El QI]] , bimaximal mixing at high energy can be 
compatible with low-energy experiments fTJ ^l or the small mass splittings can be 
generated from exactly degenerate light neutrinos [HI [TSl UHl UZl UH US- On the other 
hand, facing the high precision of future neutrino experiments, rather small RG corrections 
are important as well. For instance, deviations from ^13 = or maximal mixing 6'23 = vr/4 
are induced by RG effects |2()| El I22j also for a hierarchical spectrum. However, in most 
studies only the running of the dimension 5 operator has been considered, which is only 
appropriate for the energy range below the mass scale of the heavy singlets. 

The importance of including the effects from energy ranges above and between these 
mass thresholds when analyzing RG effects in GUT models has been pointed out in [23 
ElESlEElEZliniEllElliniEII. They are typically at least as important as the effects 
from below the thresholds since the relevant couplings, i.e. the entries of Yy, can be of 
order one, regardless of tan [3} Previous studies have investigated the RG effects above 
the see-saw scales mainly numerically. 

In this paper we derive formulae which allow to understand the running of the neutrino 
parameters above the see-saw scales analytically. We further provide a software package 
for analyzing the RG evolution (with correct treatment of non-degenerate see-saw scales) 
numerically. We apply our results to investigate consequences of the running above the 
see-saw scales for model building and leptogenesis and compare the size of RG corrections 
to the precision of future experiments. 

The paper is organized as follows: In Sec. El we review how the predictions for neutrino 
masses can be evolved from the GUT scale to the electroweak scale. Sec. El is dedicated 
to the analytic understanding of RG effects in see-saw scenarios with special emphasis on 
the range between Mqut and the highest see-saw scale. In Sec. (H we analyze the run- 
ning between the see-saw scales in more detail. Sec. El contains a brief description of the 
accompanying Mathematica packages for numerical RG analyses (a detailed documenta- 
tion is available at http://www.ph.tum.de/~rge/). In Sec. El we discuss applications to 
model building and related topics. Alternatives to the simplest see-saw scenario are briefly 
discussed in Sec. |7| Finally, Sec. El contains our conclusions. 

^Large entries of could be important in models of gauge- Yukawa unification (see, e.g., ^j), and 
may even be important for precision gauge unification in the MSSM |29j . 
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2 Running Neutrino Masses in See-Saw Scenarios 



In this section, we discuss how to obtain the RG evolution of neutrino masses, starting 
from initial conditions at a very high energy scale. ^ An important technical issue is that 
the heavy singlet neutrinos involved in the see-saw mechanism have to be integrated out 
one by one. Thus, one has to consider a series of effective theories pBl l27j. We will focus on 
the SM and the MSSM amended by three singlet neutrinos or three singlet superfields 
ivj, respectively. The discussion can be applied to other scenarios, such as multi-Higgs 
models, and a different number of singlets in a straightforward way. 
We consider the Lagrangian of the SM extended by 

= -N^YjJ^ - ^N^MN^ + h.c. , (2) 

where := (^l?^L'^l)"^ denotes the left-handed lepton doublets, (j) is the Higgs doublet 
and (p = ir'^cp* its charge conjugate. The superscript C denotes charge conjugation of 
fermion fields, and := (A^r)^. In the supersymmetric case, is replaced by the Higgs 
doublet Hu coupling to the up-type quarks. 

In order to define mass and mixing parameters as functions of the renormalization scale 
/I above the highest see-saw scale, we consider the effective light neutrino mass matrix 

mM = - Y Yjifi) M-\fi) YM , (3) 

where Y,y and M are /i-dependent. The relevant Higgs vev is f = 246 GeV in the SM and 
V = 246 GeV ■ sin jS in the MSSM.^ is the mass matrix of the three light neutrinos as 
obtained from block-diagonalizing the complete 6x6 neutrino mass matrix, following the 
standard see-saw calculation. The scale-dependent mixing parameters are obtained from 
my(/i) and the running charged lepton Yukawa matrix Ye{^). In Sec. El we are going to 
analyze the energy dependence of the parameters in the lepton sector such as neutrino 
masses, lepton mixing angles and CP phases above the highest see-saw scale analytically. 
Therefore, we will make use of the RGE for the composite quantity mj,, calculated from 
those for and M [HI 133 121 123 • It is given by 

(1 777 

IQir^ = {Ce YX + YXf + {C, Y}Y, + Y^Y,) +am, (4) 
with if: := ln(/i/ Hq), 

Ce = , a = ^ in the SM, (5a) 
C^ = C^ = l in the MSSM, (5b) 



^In the following we will refer to this high energy scale as Mqut, although it can be any other scale 
where additional new physics, apart from the heavy singlet neutrinos, has to be taken into account. 

^As indicated in Eq. Q, we do not take into account the running of the Higgs vev. In principle, v 
runs as well, so that m,y actually does not yield the physical neutrino masses. However, the evolution of v 
depends on the renormalization scheme and on the definition of the Higgs mass, see e.g. PO], so that there 
is no straightforward definition of a neutrino mass with a running vev. In any case, the mixing angles 
and phases are independent of the value of v. This definition has shown appropriate for the applications 
discussed in this paper, such as leptogenesis. 
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and (with Y^., and Y^ being the Yukawa matrices of charged leptons, down- and up-type 
quarks, respectively)'' 

«SM = -^g!-lgl + 2Tr{YX + ye^ye + ^YjY, + 3YX), (6a) 

«MSSM = -^g!~6gl + 2TT{YX + ^yIyu). (6b) 

The RGE governs only the evolution of the light neutrino mass matrix above the 
highest see-saw scale, which is given by the mass eigenvalue M3 of the heaviest singlet N^. 
For n < M3, we obtain the correct RG evolution by integrating out N^. This leads to the 
appearance of an effective neutrino mass operator 

= ^^4(f^■0)(^^0)+h.c., (7) 

where f,g & {1,2,3} are family indices and where the dot indicates the SU(2)L-invariant 
contractions. The coefficient of this operator is obtained by the (tree-level) matching 
condition^ 

'4 = 2(Fj),3M3-i(i;)3/, (8) 

which is imposed at yU = M3. This expression is specified in the mass basis for the singlets, 
i.e. in the basis where M is diagonal. Let us mention that finding the matching scale 
properly requires some care as the mass matrix M (and consequently the eigenvalue M3) 
itself is subject to the RG evolution. As a consequence, for scales below M3 the effective 
neutrino mass matrix can be described as a sum of two contributions, 

= -^{'U2YjM-'Y,) . (9) 

(3) 

The 2x3 Yukawa matrix Y^ is obtained by simply removing the last row of Y^ in the 

(3) 

basis where M is diagonal. The 2x2 mass matrix M is found from M by removing the 
last row and column. By construction, is a continuous function of the renormalization 
scale. The RG evolution of the second term on the right-hand side of Eq. is governed 

(3) 

by Eq. with Y^, replaced by Y^. The running of the first term, on the other hand, is 
determined by the RGE gZI 

1 (3) 

= {C,YX + CuYfYy''^ + ''^iCeyX + CuWt)+a^^^^ (10) 
with Ce and as in Eqs. © [SI ESI EEl EZI , and 

«SM = ~3gl + 2TT(Yft + YX + ^YjY, + 3YX^)+X, (11a) 

ttMSSM = - 6gl + 2 Tr + 3FJFJ . (lib) 
5 



^We use GUT charge normalization for the gauge coupUng gi. 

^We do not discuss finite threshold corrections, which arise due to the fact that the singlet neutrinos 
do not decouple abruptly |33j . The resulting uncertainty in the low-energy results is typically not larger 
than that due to two-loop effects. In the REAP software package described in Sec. [SJ the corrections can 
be implemented approximately by integrating out N^^ slightly below M3. 
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Figure 1: Validity ranges of the effective theories (EFTs) in the renormahzation scale ^. At a scale close to 
the mass thresholds Mi, the EFTs are related by matching conditions. Although we show this illustration 
for 3 heavy singlets, it is straightforward to generalize it to an arbitrary number (cf. |27j). 



One can now evolve the effective neutrino mass matrix down to the scale M2 and 
repeat the matching procedure there. From integrating out A^^ at /i = M2, the Yukawa 
matrix gets further reduced and the effective neutrino mass operator receives an additional 
contribution. After a subsequent RG evolution to /i = Mi, the procedure is repeated for 
A^^. The emerging effective theories, as well as the quantities relevant to neutrino masses 
in each of them, are illustrated in Fig. ^ 

In summary, the running of the effective neutrino mass matrix uiv above and between 
the see-saw scales is given by the running of two parts. 




where n labels the effective theory (cf. Fig. Q). In the SM and the MSSM, the 1-loop 
/3-functions for my in the various effective theories can be summarized as 

(n) 

„ dA^ / i ("),(") \ T ("•) (") / i (").("•) N (n) (") 

16772 — = [CeYX + CyYX) X + x{CeYjY, + CyYX)+OixX , (13) 

(") („) (") ("•) ("■) 

where X stands for k or for 2Yj; M Y^, respectively. The coefficients Cj and are listed 

in Tab. H 
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Table 1: Coefficients of tlie /3- functions of Eq. (|13() . which govern the running of the effective neutrino 
mass matrix in minimal see-saw models. 



3 Analytic Understanding of the RG Evolution 

The methods of [311^11131201 can be used to derive differential equations for the running of 
the neutrino masses, mixing angles and CP phases in the see-saw scenario. In this section, 
we concentrate on the full theory above the highest see-saw scale. The corresponding 
differential equations for the running below the see-saw scales have been discussed in 
jlOl EHl 120] • We abbreviate the flavour- dependent terms in the RGE (jH) by 

P:=CeY^Y, + C,YjY,. (14) 

Due to the appearance of the neutrino Yukawa couplings, the running depends on more 
parameters than below the see-saw scale. In particular, since the see-saw formula does 
not allow to determine Yi, uniquely from the light neutrino mass matrix, the running is no 
longer determined by (the RG extrapolation of) low-energy parameters only. Moreover, 
Y^Ye and YjY^ are not simultaneously diagonalizable in general. As a consequence, the RG 
evolution generates off-diagonal entries in the charged lepton Yukawa couplings, even if one 
starts in a basis where they are diagonal (cf. the RGEs in App. This is also different 
from the situation below the see-saw scale and makes the results more complicated. 

In a given basis, Y^Ye and 171,^ can be diagonalized by unitary matrices, Ue and U^, 
respectively. The lepton mixing matrix is given by f/jviNS = U^U^- Keeping the basis flxed, 
both matrices change with the renormalization scale, so that the RGEs of the mixing 
parameters consist of two parts, one coming from the RG change of Ue, and the other 
from the change of f/j,. We will refer to these as Ug and Ui, contribution in the following.® 
Further details and the derivation of the formulae are given in App. El 

We will flrst discuss the U,y contribution, which is often dominant. An important result 
is that in the RGEs above the see-saw scale, the same mass squared differences appear in 
the denominators as below the see-saw scale, so that 

A^i3,Ae23 oc -r^, (15b) 



^One might wonder whether it is possible to simplify the situation by working in the basis where P 
is diagonal. This is not the case, since the Ue contribution depends on a different linear combination of 
YjY, and Kjv . 
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where, as usual, Aml^^ := ml — ml and Am^^^j := m^ — m\J Thus, 612 and the phases 
generically still run faster than 6*13 and 6*23. Besides, the running is suppressed by a strong 
normal mass hierarchy, as it is the case below Mi. For the unphysical phases*, we find a 
generically larger change A5e oc l/Amg^i, while A5^, A^,- oc l/Am^^^^. 

Often, the evolution will be dominated by a single element of P. Then, the derivatives 
of the masses and mixing parameters are given by this element times the corresponding 
entry in the tables of Sec. 13.31 and App. O We will discuss an example in Sec. 16.11 Of 
course, if several entries of Pfg are relevant, one obtains the analytic description by simply 
adding up their contributions. The tables are given in the basis where Ye is diagonal and 
where the unphysical phases in the MNS matrix are zero (cf. Apps. IB. H and IB.5|) . In order 
to keep the expressions reasonably short, we only present the first order of the expansion 
in the small CHOOZ angle 6'i3. We furthermore use the abbreviation 

C := ^ . (16) 

Its current best-fit value is C ~ 0.038 [H]. Note that this value is measured at low energy. 
It can change significantly, if the running of the mass eigenvalues is not a simple rescaling. 

The tables in the appendix show that the numerators of the RGEs are of the order 
of mf in the generic case, i.e. if there are no significant cancellations. Then, the generic 
enhancement and suppression factors given in Tab. El yield a first estimate of the RG 
change of the mixing angles. In particular, they allow to understand analytically when the 
evolution is enhanced or suppressed compared to the naive estimate 

where Pfg is assumed to dominate the running and is the corresponding see-saw scale. 
The analogous factors for the CP phases are given in Tab. El The size of quantum correc- 
tions can thus be estimated by multiplying A9f^"^ with the corresponding enhancement 
or suppression factor. As the mass hierarchy is weaker in the neutrino sector than in the 
quark sector, the change of the mixing parameters in the MNS matrix is larger than that 
of the ones in the CKM matrix. 

The RG evolution can deviate significantly from the generic estimate, if cancellations 
occur. For example, for non-zero (^1 — 932, the running of 6*12 usually gets damped (as it is 
the case below the see-saw scales |12|)- Such effects can be understood from the complete 
formulae in App. O However, care should be taken when estimating the RG effects for 
special phase configurations with extreme cancellations, such as — ^2 = ^r, as terms 
proportional to ^13 (which are neglected in our formulae) can become important then. 



3.1 Running of the Mixing Angles 

From the generic enhancement and suppression factors for the evolution of the solar angle 
in Tab. |H we see that all terms in 612 are enlarged by m^ / Am^^^ for quasi-degenerate 

^For specific textures, this observation has been made in ^JIHl- The result can also be obtained by 
using the formulae of "39'. 

^The term "unphysical phases" is somewhat misleading here, since the distinction between physical 
and unphysical parameters is not completely trivial in the full theory, cf. App. IB. 51 
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Table 2: Generic enhancement and suppression factors for the evolution of the angles, yielding an estimate 
of the size of the RG effect. The table entries correspond to the terms in the mixing parameter RGEs 
with the coefficient given by the first column. A '1' indicates that there is no generic enhancement or 
suppression, 'd.' stands for a degenerate neutrino mass spectrum, i.e. ISjv?^^^ <C ^ rn^ ~ mg m? . 
'n.h.' denotes a normally hierarchical spectrum, i.e. mi <C m.2 ^ mg, and 'i.h.' means an inverted 
hierarchy, i.e. 7713 <C mi < m2. 

masses. Thus, there will be large RG effects, if the different terms do not cancel each 
other. The term involving Im P32 is an exception, because its leading order is proportional 
to ^13, so that it only plays a role in special cases. In the case of a strong normal hierarchy, 
there is no enhancement. However, for a moderate hierarchy where the square of the 
lightest neutrino mass is small compared to Am^^^^ but larger than Am^Qj the running is 
still enhanced by m\/ Aml^i. This is similar for an inverted hierarchy, where the evolution is 
generically enhanced by C^^, because the masses mi and m2 are almost degenerate. Thus, 
the RG change of 9i2 is generically large for an inverted hierarchy and for a degenerate 
spectrum, and small for a normal hierarchy. This conclusion is unchanged compared to 
the region below the see-saw scale. 

The enhancement and suppression factors of 6'i3 are similar to those of ^23- The evo- 
lution of both angles does not depend on Pu for ^^13 = 0. The terms proportional to the 
other Pfg are enhanced by m^/Am^|.j^ in the degenerate case, so that we expect significant 
effects here as well. However, as already mentioned, they are usually smaller than those for 
9i2- For both hierarchical spectra, the running is slow. For a diagonal P and an inverted 
hierarchy with = 0, ^13 does not run at all, if it vanishes at some energy, as it is the 
case below the see-saw scale jlHj. However, this is no longer true if P21 or P31 is non-zero. 

As far as the dependence of the RGEs on the mixing parameters is concerned, we find 
from Tab. 1121 that the terms in the RGEs which are proportional to the diagonal elements 
of P exhibit basically the same behavior as the RGEs below the see-saw scale |20j. The 
running of 612 and 6'23 is damped by non-zero Majorana phases, while the situation is more 
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complicated for ^13. In particular, the value of the Dirac phase in the case ^13 = is 
determined by the condition that 6 remain finite. Additionally, the running is suppressed 
if the mixing angles are small, as it is the case in the quark sector. (This is another reason 
why the leptonic mixings run faster than the quark mixings 

If the diagonal elements are equal, their contributions to the RGEs cancel exactly. 
This follows from the fact that the mixing angles do not change under the RG, if P is 
the identity matrix and thus does not distinguish between the flavours. Of course, this 
statement holds also for the RGEs of the CP phases. It provides a consistency check for 
the results. 

Interesting new effects occur for non-zero off-diagonal elements in P. Some of their 
coefficients in the RGEs do not vanish for vanishing mixings, so that non-zero mixing angles 
are generated radiatively. Because of this, it is possible to reach low-energy parameter 
regions that are compatible with experiment even if the neutrino mass matrix is diagonal 
at the GUT scale JO]- This is in striking contrast to the region below the see-saw scale 
and to the quark sector. The terms proportional to the real parts of the off-diagonal Pfg 
exhibit the same dependence on the Majorana phases as the diagonal elements. Some of 
them are suppressed for large angles 612 and 623- For example, the ReP23 contribution 
to 623 vanishes for maximal atmospheric mixing. The influence of the imaginary parts 
has quite a different dependence on the mixing parameters, in particular on the Majorana 
phases. The corresponding terms become maximal for non-vanishing phases, for instance 
for — (y92 = 7r/2 in the case of 612. Thus, the usual damping of the running by non- 
zero Majorana phases does not always take place above the see-saw scales. However, the 
maximal damping for ipi — ip2 = tt (or 99^ = vr in the case of ^23) still occurs, since the 
coefficients of Im Pjg are zero then. Some examples for the running with large imaginary 
entries in P will be discussed in Sec. 16.41 

3.2 Running of the Phases 

The CP phases show a fast running in general. The corresponding generic enhancement 
and suppression factors are given in Tab. El As for the RGE of the Dirac phase S, there is 
always a term proportional to ^^3^, which is further enhanced for a degenerate spectrum. 
This implies that the running of 6 is in general significant for small ^13, irrespectively of 
the hierarchy.^ For 613 = 0, 6 and 6 are undefined. However, it is possible to define an 
analytic continuation yielding a smooth evolution ,2^. In addition, for the degenerate or 
inversely hierarchical spectrum, the running of 6 gets enhanced by terms proportional to 
m^/Amg^j or respectively. The coefficients of Pfg in 6 are given in Tab. 1131 from where 
one obtains the RGE as S = e^^^6^-^'> + 5(0) + ^(^13). 

The situation is similar for the Majorana phases. By the same reasoning as for the 
running of the solar angle, the generic RG effects are large for degenerate masses and for an 
inverted hierarchy, while they are suppressed for a strong normal hierarchy. The coefficients 
of Pfg in (fi are given in Tab. El These formulae are also important to understand the 
evolution of the mixing angles in some cases. An example will be discussed in Sec. 16.41 

The evolution of the Majorana phase difference is governed by a simple equation, 

^Note, however, that in measurable quantities S appears always in combination with sin0i3, so that 
the RG change of predictions for experiments may not be significant. 
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Table 3: Generic enhancement and suppression factors for the evolution of the CP phases, yielding an 
estimate of the size of the RG effect. The table entries correspond to the terms in the mixing parameter 
RGEs with the coefhcient given by the first column. A '1' indicates that there is no generic enhancement 
or suppression, 'd.' denotes a degenerate neutrino mass spectrum, i.e. lS.rn^^^ <C rn\ ^ m| ~ m§ ~ ni?. 
'n.h.' denotes a normally hierarchical mass spectrum, i.e. mi ^ m2 ma, and 'i.h.' means an inverted 
hierarchy, i.e. <^ mi < 7712- 
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Table 4: Coefficients of Pfg in the slope of the Majorana phase difference for 6*13 = 0. The abbreviations 
Sij and Qfj depend on the mass eigenvalues and phases only, and enhance the running for a degenerate 
mass spectrum since they are of the form /^(mi, m^, (pi, (p2)/(jnj — mf). They are listed in Tab. ^2 We 
use the abbreviations Cy = cos 6ij and Sy = sin 9ij (cf. Ann.rOTl. 
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Table 5: Coefficients of Pfg in the slope of the mass eigenvalues for 6113 — 0. 

which can be read off from Tab. El It indicates strong running, since the slope is still 
inversely proportional to Amg^j. However, in the case of equal Majorana phases, only the 
imaginary entries in P and terms proportional to 6'i3 contribute to the running. Besides, 
the contribution proportional to the real parts is suppressed for large solar mixing. 

If Y^Yi, is close to the identity matrix, its contribution to the running is very small, 
since the terms proportional to the diagonal entries cancel approximately. Then, only the 
contribution from Y^Yf. remains, so that the evolution above the see-saw scales is essentially 
the same as below. However, many GUT models suggest a hierarchical structure for Fj, 
like for the other Yukawa matrices. Then the main contribution will be due to P33 and the 
next-to-leading contribution will be from ReP32, if Y^Yy is almost diagonal in the basis 
with diagonal Y^Ye. Thus, the phase difference tends to decrease while running down,^° 
as it is the case below the see-saw scales. 

3.3 Running of the Masses 

Below the see-saw scales, the evolution of the mass eigenvalues is, to a good approximation, 
described by a universal scaling caused by the flavour-independent part of the RGE [301 
123 120] • This flavour-independent term, however, becomes smaller at high energies. In the 
MSSM, it can even cross zero at some intermediate scale. Therefore, the flavour-dependent 
terms play a more important role above the see-saw scales, the more so they can be larger 
if the entries of Yy are order one. 

We list the coefficients in the slope of the mass eigenvalues and of the Am^s in Tab. El 
and Tab. IHl respectively. Clearly, the RGE for each mass eigenvalue is proportional to the 
mass eigenvalue itself. As a consequence, the mass eigenvalues can never run from a finite 
value to zero or vice versa. In other words, the rank of the effective neutrino mass matrix 
is conserved under the renormalization group. In contrast, the mass squared differences 
can, in principle, run through zero. This, however, requires a very high value of mi. 

The flavour-independent term in the MSSM is subject to large cancellations (cf. Eq. 
(I6b|) ). Note that the running of the mass eigenvalues strongly depends on the top Yukawa 

^°More accurately, it runs away from vr and towards either or 27r, i.e. — decreases for — 1^2! < 
and increases for — (^2! > tt. 
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Table 6: Coefficients of Pfg in the slope of the mass squared differences for 613 = 0. 

coupling yt, since the term a contains 6y^, and on the gauge couphngs, which run differently 
for different SUSY breaking scales. This could, at least partially, explain why there exist 
mutually inconsistent numerical results for the scaling of the mass eigenvalues below the 
see-saw scales 1201 ESI EE] • 

Between and above the see-saw scales, the running is strongly influenced by the neutrino 
Yukawa couplings. In particular, depending on the size of the entries, omssm can 
turn negative or not. For order one entries, it typically stays positive. However, in 
such a situation, q;mssm becomes small so that P can dominate the running. Consider, 
for instance, the case of a dominant P33 entry. Here, the coefficient of m2 is enhanced 
compared to the rhi coefficient by (m2/mi) cot^ 6*12 (cf. Tab. El). In many cases 6*12 is at 
high scales much smaller than its low-energy value, so that 1712 runs much faster than mi. 
As a consequence, Amg^j can be significantly enhanced even for not too degenerate spectra. 
A relatively drastic example is shown in Fig. |21 Clearly, the discrepancy in the scaling of 
Am^Qj and Aml^^ stems from the flavour- dependent terms P. As tan/5 is large in this 
example, the P33 induced terms cause important effects already below the see-saw scale. 
The dominant effect, however, is the running in the range M3 < fi < Mgut, i-e. over less 
than two orders of magnitude. By inspecting the tables, we find that analogous features 
are present if other elements of P are large. In particular, one can enhance the evolution 
of 

A?^atm well. Therefore we expect many models which predict realistic values for the 
masses at tree level to be ruled out by several standard deviations due to KG effects. 

If, on the other hand, the eigenvalues of YjY^, are much smaller than 1, omssm typically 
flips its sign. The entries of P are now small if tan (3 is small, and for large tan (3 they 
are dominated by YjYe. Hence, for small tan /5, omssm still dominates the running of the 
masses (away from its zero point). In contrast, for large tan f3, the contribution of P (being 
now dominated by Y^Yg) is of similar importance, as it is the case for the running of the 
effective neutrino mass operator k at high energies. Since a can be negative at scales close 
to the GUT scale now, the contributions from the diagonal entries in P can decrease the 
RG effects. The off-diagonal entries again can both increase and decrease them. 

Finally, let us mention that since the terms in involving the imaginary part of P 
are proportional to sin 6*13, they do not contribute in the approximation of vanishing 6*13. 
Clearly, in the SM, a dominates the running if Y^, is small. 
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Figure 2: Example where the flavour-dependent terms dominate the running of the mass eigenvalues for 
M-i<[i< A/gut in the MSSM. We use Yy = diag(0.02, 0.1, 1) and mi = 0.04 eV at the GUT scale as well 
as a SUSY breaking scale of 200 GeV and tan /? = 50. M is chosen such that the low-energy parameters are 
compatible with experiment. The different gray-shaded areas indicate the ranges of the effective theories 
(cf. Fig.UJ. 



3.4 Ue Contribution to the Running 

As mentioned in the beginning of this section, the RGE for Yg contains non-diagonal 
terms above and between the thresholds, so that there is an additional contribution to the 
running of the leptonic mixing angles and CP phases. In the see-saw scenario, the RGE 
for Yf. above M3 is given by 

Ay 

IGtt^ ^ = Ye (D, Y^Y, + D, YJY,) + Y, =:Y,F + Y, (18) 



dt 



with 



De = ^, D, = -^ in the SM, (19a) 
De = 'i, D^ = l in the MSSM. (19b) 

As usual, Ue is flavour diagonal (cf. App.EJ- The resulting contributions to the evolution 
of the angles for vanishing 61^ and <^ Vt are listed in Tab. [7| They can simply be 

added to the expressions discussed above (cf. App. 16.4]) . 

In contrast to the latter, all non-zero terms in the contribution have a generic 
enhancement factor of 1. The reason for this is the strong hierarchy among the charged 
lepton masses. As a consequence, the Ue contribution is negligible compared to the U„ 
contribution, if the relevant factor in Tab. El is much larger than 1. If it is close to 1, 
both contributions are generically of the same order of magnitude. The f/g contribution 
can even be dominant if the factor is small. This is also possible, if cancellations occur 
between the leading-order terms in the RGEs. 

To get a feeling for the size of the effects discussed in this section, let us consider a 
rough estimate. We assume that the running is linear on a logarithmic scale, that it is 
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Table 7: Coefficients of Ffg in the Ue contribution to the slope of the mixing angles for O13 — and 

dominated by a single entry y in Y^, which is related to the light neutrino mass rris and 
the see-saw scale M3 by = yms' ^'^^^ relevant term in Tab. [3 is of the order 
of 1. Then we find 

lA^^I ~ l^^l Duv' (0.O27 + O.OO6 ln!!^^ZH^) . (20) 

Thus, the change is small, but it can still be relevant in the context of precision studies 
(e.g. the change of 613), if y is large. 

4 Running between the See-Saw Scales 

Between the see-saw scales, the singlets are partly integrated out, which implies that only 
a (n — 1) X 3 submatrix of the neutrino Yukawa matrix remains. Therefore, we expect 
that the running between the thresholds caused by the neutrino Yukawa matrix can differ 
significantly from the running above or below them. 

We now discuss the running due to the terms in the /3-functions with a flavour structure 
proportional to the unit matrix. Below and above the see-saw scales, they only cause a 
common scaling of the elements of the neutrino mass matrix and thus leave the mixing 
angles and phases unchanged. Between the thresholds, however, the effective neutrino 

mass matrix consists of the two parts k and 2Yj M Yy, as shown in Eq. (fT^ . Here, the 
mixing angles and phases change in general, unless both parts are scaled equally. From 

(n) (n) 

tablelH we see that in the SM, the /3-functions and (32yjm-'^y^^ have different coefficients 
in the terms proportional to the gauge couplings and to the Higgs self-couphng [2Z1- This 
difference can be understood by looking at the corresponding diagrams of the "full" and 
the effective theory. For instance, the diagram for the correction to the effective vertex 
proportional to A and its counterpart with the heavy singlet running in the loop are shown 
in figure El Diagram (a) is UV divergent, whereas diagram (b) is UV finite. We thus get no 
contribution proportional to A for the /3-function of the composite operator. The situation 
is similar for some of the diagrams corresponding to the vertex corrections proportional to 

the gauge couplings. Thus, in the SM, the KG scaling of the two parts k and 2 Y^ M Y^ 



13 




Figure 3: Figure (a) shows the diagram which gives the contribution proportional to the Higgs self- 
coupUng in the /3-function of the neutrino mass operator in the SM. Figure (b) shows its finite counterpart 
with the heavy singlet running in the loop. The gray box labeled by k* corresponds to the contribution 
to the effective neutrino mass operator from integrating out the heavy singlet . 

of the effective mass matrix between the thresholds, caused by the interactions with trivial 
flavour structure, is different. This implies a running of the mixing angles and CP phases in 
addition to the running of the mass eigenvalues.^^ This effect can even give the dominant 
contribution to the running of the mixing angles, as for instance in the example shown in 
figure m (from ) . 

(n) (n) 

Due to the non-renormalization theorem in supersymmetric theories, (3^ and (32YJ'm-^y^ 
are identical in the MSSM (see Tab. [T]on p. EJ, so that we can use the RGEs of Sec. El 
between the see-saw scales as well. In particular, the enhanced running between the 
thresholds due to terms with a trivial flavour structure does not occur. Of course, the 
heavy degrees of freedom have to be integrated out first, i.e. all parameters have to be 
replaced by the effective ones between the thresholds. 

5 Mathematica Packages for Numerical RG Analyses 
5.1 Numerical Solution of the RGEs 

The Mathematica package REAP (Renormalization Group Evolution of Angles and Phases) 
numerically solves the RGEs of the quantities relevant for neutrino masses, for example 
the dimension 5 neutrino mass operator, the Yukawa matrices and the gauge couplings. 
The /3-functions for the SM, the MSSM and two Higgs doublet models with Z2 symmetry 
for FCNC suppression (2HDM) with and without right-handed neutrinos are implemented. 
In addition, the same models are available for Dirac neutrinos. New models can be added 
by the user. The heavy singlet neutrinos can be integrated out automatically at the 
correct mass thresholds, as described in SecEl^^ The software can also be applied to type 
II see-saw models as long as one only considers the energy region below the additional 
see-saw scale Ma, where the new physics such as Higgs triplets only leads to another 

"To see this, let us assume that [/^('k + 2YJM^'^Y^) U is diagonal. Then U'^{aK + b2YjM-'^Y^) U 
is in general only diagonal if a = 6 (common scaling). 

^^We do not consider SUSY threshold corrections they are usually numerically less important 

EH!- 
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Figure 4: Running from maximal solar mixing at A/qut to the experimentally preferred angle of the 
LMA solution. The figure shows an example in the SM with a negative CP parity for m2 and a Yukawa 
matrix = 0.5 • diag(£^, e, 1) at Mgut with e = 3.5 • 10~^ and a normal mass hierarchy (from 55)- The 
lightest neutrino has a mass of 0.004 eV (at low energy). The gray-shaded areas illustrate the validity 
ranges of the effective theories emerging from integrating out the heavy singlet neutrinos. 



contribution to the effective neutrino mass operator. Tlie pack;age can be downloaded 
from http : / / www . ph . turn . de/~rge/REAP/ , Matliematica 5 is required. 



5.2 Extraction of Mixing Parameters from Mass Matrices 

Tlie package MixingParameterTools" ("MPT) allows to extract the physical lepton masses, 
mixing angles and CP phases from the mass matrices of the neutrinos and the charged 
leptons. Thus, the running of the neutrino mass matrix calculated bv IREAPl can be trans- 
lated into the running of the mixing parameters and the mass eigenvalues. For the def- 
inition of the mixing parameters, see App. lA.ll and the documentation of the package. 
|HixingParameterTools can also be useful as a stand-alone application in order to study 
textures without running, and it is not bound to the analysis of neutrino masses only 
but may be used for quark and superpartner mass matrices as well. Therefore, it can be 
obtained separately from IREAPl at http : //www . ph . turn . de7~rg e/MPT/1 



5.3 Example Calculation 

The following simple example demonstrates how to use the Mathematica packages to cal- 
culate the RG evolution of the neutrino mass matrix in the MSSM extended by three 
heavy singlet neutrinos. Of course, further documentation is provided together with the 
packages. 

1. The package corresponding to the model at the highest energy has to be loaded. 
All other packages needed in the course of the calculation are loaded automatically. 
(Note that ' is the backquote, which is used in opening quotation marks, for example.) 

Needs [ "REAP ' RGEMSSM ] 
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2. Next, we specify that we would like to use the MSSM with singlet neutrinos and 
tan/9 = 50. Furthermore, we set the SUSY breaking scale to 200 GeV and use the 
SM as an effective theory below this scale. 

RGEAdd["MSSM",RGEtan\ [Beta]->50] 
RGEAdd["SM" .RGECutof f->200] 

3. Now we have to provide the initial values. For instance, let us set the GUT-scale 
value of 9i2 to 45° and that of the first Majorana phase to 50°. Besides, we use a 
simple diagonal pattern for the neutrino Yukawa matrix and the default values of 
the package for the remaining parameters. 

RGESetlnitial [2*10" 16 , 

RGE\ [Theta] 12->45 Degree , RGE\ [Phi] l->50 Degree , 
RGEY\ [Nu]->{{1,0,0},{0,0.5,0},{0,0,0.1}}] 

4. RGESolve [low, high] solves the RGEs between the energy scales low and high. The 
heavy singlets are integrated out automatically at their mass thresholds. 

RGESolve [100,2*10^16] 

5. Using RGEGetSolution [scale, guantity] we can query the value of the quantity 
given in the second argument at the energy given in the first one. For example, this 
returns the mass matrix of the light neutrinos at 100 GeV: 

MatrixForm [RGEGetSolution [100 , RGEM\ [Nu] ] ] 

6. To find the leptonic mass parameters, we use the function MNSParameters [m,^ jFe] 
(which also needs the Yukawa matrix of the charged leptons). The results are given 
in the order {{Ou, 6*13, 623, S, Se, S^, Sr, (pi, (P2}, {mi, m2, ma}, {ye, y^, yr}}- 

MNSParameters [ 

RGEGetSolution [100, RGEM\ [Nu]] , RGEGetSolution [100, RGEYe]] 

7. Finally, we can plot the running of the mixing angles: 

Needs ["Graphics 'Graphics' "] 

mNu [x_] : =RGEGetSolut ion [x , RGEM\ [Nu] ] 

Ye[x_] :=RGEGetSolution[x, RGEYe] 

\ [Theta] 12 [x_] : =MNSParameters [mNu [x] , Ye [x] ] [ [1 , 1] ] 
\ [Theta] 13 [x_] : =MNSParameters [mNu [x] , Ye [x] ] [ [1 , 2] ] 
\ [Theta] 23 [x_] : =MNSParameters [mNu [x] , Ye [x] ] [ [1 , 3] ] 
LogLinearPlot [{\ [Theta] 12 [x] , \ [Theta] 13 [x] , \ [Theta] 23 [x] } , 
{x,100,2*10~16}] 
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6 Applications 



We now apply the analytical and numerical tools described in the previous sections to 
some specific cases with interesting RG effects above, between and below the see-saw 
scales within the conventional see-saw scenario. 



6.1 RG Effects for a Dominant (1^)33 



Many unified models relate the Yukawa couplings of the different charged fermions and 
the neutrinos, e.g. ~ or Y^, ~ Yg. For the charged fermions, the quantities accessible 
through observation are Y'^Y, where Y denotes the corresponding Yukawa matrix. It is 
convenient to work in the basis where Yj^Yu and Y^Ye are diagonal and positive, and the 
diagonal entries are ordered ascendingly. In this basis, all three combinations Y'^Y have 
a dominant 33 entry. In this subsection, we shall assume a similar pattern for YjY^, 
i.e. {YJYi,)33 ^ t/g ^ (YjYi,)ijji33. Given such a hierarchy for YJY^, the RG corrections 
A9i3 := 6'i3(Msusy) — 6'i3(Mgut) and A923 can be approximated by 
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where denotes the mass scale of the heavy neutrino (s) with the large Yukawa cou- 
plings.^^ To obtain these results, we read off the RGEs from Tab. fT^ and integrated them 
with the approximation of constant coefficients. This is reasonably accurate, since the 
running of 6'i3 and 6'23 is almost linear on logarithmic scales [20] -^"^ 

In the SM, the term proportional to y1 is negligible, since the Yukawa coupling is not 
enhanced by tan (3. However, the y1 contribution can be large, and it is not suppressed for 
small tan j3. Furthermore, except for and M^,, only (the RG extrapolation of) low-energy 
parameters enter the expressions (PT|) and ()22|) . 

In the case of the solar angle, the running is strongly non-linear when the RG change 
is large. Then, the approximation used in the above equations does not yield reliable 
results. Even by integrating the RGE (assuming 9i2 to vary but the other parameters to be 
constant), one arrives at an expression which does not represent an accurate approximation 
in many cases because of the running of Am^^^j. Nevertheless, an inspection of the RGE 
reveals several qualitative features of the running such as the damping infiuence of the 
phases, as discussed in Sec. 13.11 

The running of the Majorana phases may be regarded as encouraging for the prospects 
of neutrinoless double {3 decay experiments: it is known that even if the mass eigenvalues 

^■^For the analytic estimates, we ignore complications due to the generically non-degenerate see-saw 
scales 27 . 

^"^A comparison with numerical calculations shows that this is unchanged in the presence of neutrino 
Yukawa couplings. 



17 



are large enough to make a discovery in future experiments possible, cancellations may 
strongly suppress the amplitude jl9]. This can directly be seen from the fact that the 
amplitude is governed by the effective neutrino mass 



/ \ I 2 2 I 2 2 iwo I 2 



2i5| 



(23) 



which is obviously suppressed if ipi — ip2 is close to tt. However, for dominant P33, the 
difference of Majorana phases is driven away from tt at low energies due to RG effects (cf. 
the discussion in Sec. I3.2|) . This implies that cancellations tend to be avoided. Note that 
the contribution from YjYf., which persists below the see-saw scales, increases the effect 

EI 



6.2 Neutrino Yukawa Couplings with Two Large Entries 

As another example, let us assume that the neutrino Yukawa matrix contains two domi- 
nant entries, {Y,y)^s ^ G~^''{Y,y)s2 ~ with an arbitrary phase 7, as it is the case in many 
models where the large atmospheric mixing angle emerges from in the basis where Y^ 
is diagonal. Then {YjY^)^^ ^ {YjY^)22, which causes a cancellation between the contribu- 
tions proportional to these terms in the RGEs of 6*13 and 6'23. Thus, using the same linear 
approximation as in Sec. 16.11 we obtain the changes 
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The change proportional to the real part of P32 vanishes for maximal atmospheric mixing. 
Hence, the neutrino Yukawa couplings only contribute significantly to the running of ^13 
in this case, if {Yy)s2 has a large imaginary part and if the CP phases are not close to or 
TT. In A^235 they always play a role by inducing off-diagonal elements in Y^Yf., which leads 
to the last term in Eq. . This term is actually dominant in the case of CP conservation 
and small tan /3. 
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6.3 RG Corrections and Precision Measurements 

In this section, we will estimate the order of magnitude of RG effects in see-saw models 
and compare it to the precision of future measurements of neutrino mixing (see also [211 
HU EHl m] for related works). We shall first consider the effects of a large P33 as an 
example. For instance, P33 can be generated from the entry (y^)33. Note that this is only 
an example. RG effects from different structures of Y^, can be understood and estimated 
using the analytic formulae of Sec. El Graphically, the RG corrections caused by P33 in 
the MSSM with tan (3 = 20 are illustrated in Fig. El We have assumed the initial values 
6*13 = 0, 6*23 = vr/4 and 9i2 + 9c = 7i/4: (where 9c is the Cabibbo angle) at high energy, 
which may be especially interesting from a theoretical point of view |221 EHl E^- The 
changes of ^13 and ^23 have been calculated from the approximations ()21|) and ()22|). We 
would like to stress that the mass squared differences are running quantities as well and 
taking them as constant, as it was done in Eqs. (PT|) and ()22|1 . restricts the accuracy of 
the estimates. For producing the plots in Fig. El we have used the values of Am^^m and 
Am-gQi at = 10^^ GeV. For the considered parameter ranges and for mt{mt) = 175 GeV 
and MsusY = 1 TeV, the mass squared differences at = 10^^ GeV are about a factor 
1.75 larger than the low-energy values. Note that their running depends sensitively on 
the value of the top mass and on the SUSY breaking scale. The change of 9i2 has also 
been determined assuming a linear running, which is possible here because only rather 
small neutrino masses and a moderate value of tan/3 are considered in the plot. We have 
used those values for the Majorana phases that do not damp the RG evolution, as well 
as best-fit values for the oscillation parameters. For the see-saw scale associated with the 
large Yukawa coupling, we have used the approximation 

M,^M33^y(n)L(m;%. (26) 

To justify this, let us reconstruct M from and mj, using the inverse of the see-saw formula 
0, M = —^Yym~^Yj, for a dominant entry {Y^)^^ in Yy and not too large neutrino 
masses, mi < 0.1 eV. In this case, one can see from = Uy diag(m]~^, \ WI3 ^) f/J 
that all entries of the inverse light neutrino mass matrix are usually of the same order 
of magnitude. Consequently, M33 is dominated by the term proportional to (yj/)|3, i.e. 
the one given in Eq. ()26|) . Furthermore, M33 is the dominant entry in M, so that it is 
approximately equal to the largest eigenvalue M3 = M*. 

We find that the RG changes are comparable to the sensitivities of planned precision 
experiments (cf. Tabs. El and EI) in the shaded parts of the parameter space, providing a 
reason to be optimistic about the potential of these experiments to find interesting results 
and to constrain model parameters. Compared to the change due to the charged lepton 
Yukawa couplings alone pUj, the gray-shaded regions are expanded, since the contribution 
from the neutrino Yukawa couplings has the same sign in the case we considered. For a 
very strong mass hierarchy, we find very small RG effects in our example. One reason for 
this is the decrease of the enhancement factors in the RGEs, as discussed in Sec. but 
this is not the main effect. What is more important is the increase of M^,. From Eq. ()26|) we 
find that it is roughly proportional to m^^ for a strong hierarchy, so that it becomes close 

^^Only for a narrow range in rni and a large difference of the Majorana phases, a suppression of the 
element (m^^)33 is possible. Then, Eq. 1)26(1 may not be a good approximation. 
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Figure 5: Estimated RG corrections to 613 = 0, ^23 = 7r/4 and 612 + Oq — vr/4 with a large P33 in the 
MSSM with tan/3 = 20, AfsusY = 1 TeV and a normal neutrino mass ordering. For instance, P33 can be 
generated from the entry (Y^)^^ in the neutrino Yukawa matrix, which was assumed here. The running 
between the electroweak and the GUT scale has been calculated using the approximate formulae (|21|1 and 
(|22|l . For producing the plots we have used Am^tm and Am^Qj a.t fj, — 10^^ GeV, which, for the considered 
parameter ranges, are about a factor 1.75 larger than the low energy values. In Fig. (a) and (c) the CP 
phases have been set to zero, and in Fig. (b) ipi — and (p2 = was assumed, leading to un-suppressed 
running. Besides, the initial condition ^13 = as well as the best-fit values for the remaining parameters 
have been used. 
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Current 


Beams 


D-CHOOZ 


T2K+NuMI 


Reactor-II 


JPARC-HK 


NuFact-II 


0.14 


0.061 


0.032 


0.023 


0.014 


10-3 


6 X 10-5 



Table 8: Current and expected sensitivities for sin^ 2^13 at the 90% CL |55l I5f)l 1^ . The entry "Beams" 
includes the conventional beam experiments MINOS, ICARUS and OPERA. The last entry refers to an 
advanced stage neutrino factory with experiments at two different baselines. The sensitivity of a first 
stage neutrino factory ("NuFact-I") is similar to that of JPARC-HK. For a description of the experiments 
and the assumptions used in the analysis, see [551 1561 15 7| and references therein. The numbers should 
be treated with some care, since they depend on the true values of the other oscillation parameters, in 
particular Am^^j^^. 



Current 


Beams 


T2K+NuMI 


JPARC-HK 


NuFact-II 


0.16 


0.1 


0.050 


0.020 


0.055 



Table 9: Current and expected sensitivities for |0.5 — sin^ ^23! |22j . The numbers are the minimal values 
required to exclude maximal mixing at the 90% CL. "Current" is the current limit from SuperKamiokande 
|58j . "Beams" means conventional neutrino beams. See ^^'^ references therein for a description of the 
experiments and the analysis methods. As in Tab. [SJ the results depend on the true values of the other 
oscillation parameters. 



to or even larger than Mqut- Consequently, the RG effects from (Y^)^^ become negligible, 
and we are left with the change proportional to y'^. This change is small here, since we are 
using a moderate value of tan/3 = 20. 

In order to demonstrate that RG corrections from Yi, are not necessarily negligible for 
a strongly hierarchical spectrum, let us consider another example, where two elements of 
Yj, are large. The evolution of the atmospheric mixing angle and mass squared difference 
is shown in Fig. IHlfor 6'23 = vr/4 at high energy in the MSSM with different values of tan f5 
and a strong normal mass hierarchy. In this example, we have taken {Y^)^^ = {Y,y)s2 = 1 
at M3 and assumed the other entries in Yi, to be small in the basis where M and Ye are 
diagonal. We have furthermore assumed that the right-handed neutrino with mass M3 
dominates in the see-saw formula, as it is the case for heavy sequential dominance (HSD) 
|59| I6()j.^^ This allows to approximately calculate M3 fa v"^ (Yy)"^^ m^^ with ma ~ \/ Am^^m 
in this case, and to consider only one see-saw scale = M3 when discussing the running. 
Eq. then simplifies to 

A^23 - ^Hiiw) V + 2VC^i2Cos^2j + ^ ln( j . (27) 

The resulting change of 6'23 is in the range of about [1°,5°]. Thus, even with a strong 
normal mass hierarchy, the change of the mixing angles can be within the sensitivity of 
future long baseline experiments. The phase yji is irrelevant due to mi = 0, and (p2 cannot 
cause a significant damping as it appears together with the rather small quantity v^- In 
Fig. |B[ it has been set to 0. 

As argued in Sec. 16.21 the running of Ug (the second term in Eq. ()27|)) cannot be ne- 
glected in this example, because the U^, contribution is strongly suppressed due to the 
cancellation between the terms proportional to P22 and P33 and the vanishing of the term 
proportional to P23 for maximal atmospheric mixing and real Y^. Even without cancella- 
tions, both contributions are generically of the same order of magnitude for hierarchical 

^^RG effects in this case have been discussed numerically in |26| . in agreement with our analytic results. 
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(a) Evolution of the atmospheric mixing angle 023 




Figure 6: Example for the running of 623, Fig. (a), and Arn^^^, Fig. (b), for a hierarchical neutrino 
spectrum. The plots show the RG evolution in the MSSM for tan/3 — 55 (solid lines), 40 (dashed lines) 
and 10 (dotted lines) with 6*23 = 45° at high energy and present best-fit values for the other parameters as 
constraints at low energy. We have used {Y^)^^ = (5^i/)32 = 1 at the see-saw scale AI3 (in the basis where 
M and Ye are diagonal) as an example (note that we use RL-convention for F^). We have furthermore 
assumed that the right-handed neutrino with mass M3 dominates in the see-saw formulae, as in heavy 
sequential dominance |59l , which allows to approximately calculate M3 from 7713 in the hierarchical 
scheme. To a good approximation, only one see-saw scale is relevant for the running in this case. The 
gray regions correspond to energies above this scale. The evolution of Ato^^j^j depends quite sensitively 
on the value of the top mass and on the SUSY breaking scale. We have used mt{mt) — 175 GeV and 
MsusY = 1 TeV. 
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Figure 7: Fast running of the solar angle despite large Majorana phases tpi = 7r/2, (/32 = in the MSSM 
with tan /3 = 30, AfsuSY — 1 TeV and a normal mass hierarchy. The evolution is dominated by the large 
imaginary part of P31, see Eq. (I28II . Further initial conditions at the GUT scale Mgut = 2 • 10^^ GeV 
were bimaximal mixing, nii = 0.08 eV, Ato^qj = 1.2 • 10^^ eV^, and Am^^j^^ ~ 4 • 10^'^ eV^. 



neutrino masses. Another lesson that can be learned from this example is that a complete 
cancellation of the running is very unlikely. Hence, we always expect RG effects to be 
comparable to the sensitivity of planned precision experiments if there are large Yukawa 
couplings and if Yi, and Ye are not simultaneously diagonal. 

6.4 Large RG Effects Despite Phases 

The main new effect above the see-saw thresholds is the appearance of off-diagonal terms in 
the Yukawa couplings. As large off-diagonal entries in the Yukawa matrices are postulated 
in a lot of fermion mass models in order to explain the large lepton mixing angles, we 
expect an important impact on the running in many cases. As mentioned in Sec. 
the effect of large imaginary entries in P is especially unusual, since their coefficients in 
the RGEs of the mixing angles 0i2 and 023 vanish for zero Majorana phases and become 
maximal if the phases or their difference equal 7r/2. Thus, a fast running is now also 
possible for large Majorana phases. A numerical example with 

/O.OOl 0\ / 0.16 0.2i\ 

>;(Mgut) = 0.01 ^ F;n(MGUT) = 0.0001 , (28) 
\-0.4i 0.5/ \-0.2i 0.25/ 

i.e. a large and purely imaginary P31 (as usual given in the basis where Y^. is diagonal and 
all unphysical phases are zero) is shown in Fig. [7| We used the MSSM with tan/? = 30, 
MsusY = ITeV, a normal hierarchy, mi = 0.08 eV, Am^^i = 1.2 ■ 10-^eV^ Am^^^^ = 
4-10-3 eV^ ipi = 7r/2, = and bimaximal mixing at the GUT scale Mgut = 2-10^^ GeV. 
Reasonable values for the low-energy oscillation parameters are reached, and Am^^j stays 
positive. The running of the solar angle from maximal mixing to smaller values is caused 
by the term proportional to Im P31 in the RGE. A negative value of Im P31 is required for 
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6i2 > (cf. Tab. IT^ . which is necessary to avoid running to the "dark side" of the solar 
oscillation parameters (corresponding to Am^^j < with our conventions). Alternatively, 
one could choose ImPsi > and exchange the initial phases, i.e. (pi = 0, ip2 = 7^/2. The 
terms proportional to the diagonal elements Pu and P33 do not play a significant role 
here, since they have opposite signs and therefore cancel approximately. The example 
demonstrates that for sufficiently large off-diagonal entries in YjY^,, it is possible to avoid 
the requirement of an inverse hierarchy of the neutrino Yukawa couplings which was found 
for diagonal YjY^ dH [121113] • 

Adding another large imaginary entry in the 32-element, 

/O.OOl 0\ / 0.16 0.2 0.2i\ 

^(Mgut) = 0.01 ^ YXiMcvT) = 0.2 0.25 0.25i , (29) 
\-0.4i -0.5i 0.5/ \-0.2i -0.25i 0.25/ 

yields a rather extreme behavior of 612, as shown in Fig. |H1 The highest see-saw scale lies 
at about 8 • 10^^ GeV here, i.e. the turnaround in the running is not a threshold effect. 
Instead, it is due to the evolution of the Majorana phases, cf. the lower plot in Fig. |H1 
Their difference initially equals tt/2 but quickly starts to increase as soon as 612 has moved 
away from 7r/4. The evolution is dominated by the term proportional to ImP^i, which 
is largest for ipi — f2 = tt. At this point, sin((^i — (^2) changes its sign, causing a sign 
change in the contributions of the imaginary parts of the off-diagonal Yukawa couplings 
to the RGE for 612. This explains the minimum in the evolution of this angle. At lower 
energies, the difference of the Majorana phases reaches a value of about 4.4 and remains 
approximately constant afterwards. From Tab. El one would expect this value to be 
closer to 27r. The difference is due to the subleading contributions to the running (the 
terms proportional to sin6'i3 and the charged lepton contribution), which become relevant 
here because of the strong damping of the leading terms. 

6.5 Leptogenesis and RG Corrections 

Leptogenesis [H^ is an attractive explanation of the observed baryon-to-photon ratio 
ub/tIj = {6.5~^q I) ■ 10~^° [02] • It typically operates at the mass scale of the lightest right- 
handed neutrino. In such a scenario, we have to deal with three scales: the GUT scale 
where the predictions for the model parameters are fixed, the scale of leptogenesis where 
the parameters have to be right for successful baryogenesis, and the low scale at which 
the parameters can be measured in experiments. In particular, one cannot use GUT scale 
parameters or experimental results directly in order to test the viability of leptogenesis in 
a given model, rather one has to take into account quantum corrections. In the energy 
range between the leptogenesis scale Mi and the electroweak scale Mew? we can consider 
the running of the effective neutrino mass operator. For relating the see-saw parameters 
at the GUT scale with the ones at Mi, the evolution above and between the see-saw scales 
has to be considered. 

"'^''This happens even if the heaviest singlet neutrino is not integrated out, i.e. even if the large Yukawa 
couplings are not removed from the theory. 
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Figure 8: Highly non-linear running of 612 and the Majorana phases in an example with large imaginary 
entries in the neutrino Yukawa matrix (see Eq. (1^ 1. We used the MSSM with tan /3 — 10, Msusy = 1 TeV 
and the following initial conditions at Mqut = 2 • 10^^ GeV: 9i2 = 023 = 7i'/4, ^13 = 0, iy9i = 7r/2, ip2 — 0, 
normal hierarchy, mi = 0.08 eV, Amf^j = 1.1 • 10^'' eV^, Am^tm = 4 • lO'^eV^ 
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6.5.1 Corrections to Decay Asymmetries and to the Neutrino Mass Bound 



The decay asymmetry for leptogenesis in the SM jnHI can be written as 




(30) 



if Ml <^ M2,M3. In the MSSM, it is a factor of 2 larger. In the case of a type II see- 
saw and for Mi <C Ma, where Ma is the mass of the SU(2)L-triplet Higgs, the decay 
asymmetry for type II leptogenesis via the lightest right-handed neutrino coincides with 
the result for the conventional see-saw [03] . In the SM or for a moderate tan f3 in the 
MSSM, the RG running from Mew to Mi leads mainly to a scaling of the neutrino mass 
matrix mj,. Including the RG effects results in an enhancement of the decay asymmetry 
for leptogenesis by roughly 20% in the MSSM and 30% - 50% in the SM [HSl I^Dl- The 
decay asymmetry can be calculated by the REAP package described in Sec. El as a function 
of energy. Thus, one can easily check if a particular high-energy model for fermion masses 
is able to produce a large enough asymmetry. Let us remark that also the running of the 
mixing angles can be very important for the calculation of the baryon asymmetry, as has 
been shown recently for non-thermal leptogenesis models f66|. 

The requirement of successful baryogenesis via thermal leptogenesis imposes constraints 
on fermion mass models and even places an upper bound on the mass of the light neutrinos 
|H7j . With respect to quantum corrections to this mass bound, it turns out that there are 
two effects operating in opposite directions, which partially cancel each other HH j IHHj: 
on the one hand, the increase of the mass scale leads to a larger decay asymmetry compared 
to the one at low energies. On the other hand, it results in a stronger washout driven by 
Yukawa couplings. Taking into account these effects and further corrections, one finds that 
the upper bound on the neutrino mass scale becomes more restrictive. 

6.5.2 Models for Resonant Leptogenesis and RG Corrections 

As an example where the running above the lowest see-saw scale can have large effects, 
we consider the RG corrections to the small mass splitting AM = |Mi — M2I for resonant 
leptogenesis [HHl EH HOI El ■ Here, the decay asymmetry is enhanced compared to Eq. (jHU)) . 
For resonance effects in the decay asymmetries to be maximal, a mass splitting of | times 
one of the decay widths (in the MSSM) 



with M2 ~ Ml := Mq, is required. Given a model for neutrino masses with such a small 
mass splitting defined at Mqut, the decay rate can be affected significantly by the RG 
evolution of the mass matrix of the heavy right-handed neutrinos from Mqut to Mi ^ M2. 
Resonant leptogenesis with exactly degenerate heavy singlets at Mqut has been discussed, 
e.g., in ^SlZnilZll- The running of M and Yi, between Mgut and Mi, taking into account 
the effects between the see-saw thresholds, can be computed conveniently using the software 
packages presented in Sec. 





22 , 



(31) 



26 



7 Alternative Scenarios 



For the examples in Sec. El we have focused on the conventional see-saw mechanism in the 
SM and in the MSSM. We now give a brief outlook on other scenarios. Some of them are 
already implemented in the software packages REAP /MPT introduced in Sec. El 



7.1 Type II See-Saw 

A generalization of the conventional see-saw is the type II see-saw [73 [701 [III 5 where 
an additional contribution to the neutrino mass matrix, e.g. from an induced vev of a 
SU(2)L-triplet Higgs, is present. Below the additional see-saw scale given by the mass 
Ma of the triplet, it can be integrated out, only leaving an additional contribution to the 
effective neutrino mass operator. The packages REAP /MPT and the analytic formulae for 
the running of the neutrino parameters can thus be applied for analyzing type II see-saw 
scenarios below M^. Above M^,, the RGEs are modified due to the additional interactions. 



7.2 Dirac Neutrinos 

At present it is not known whether the nature of neutrino masses is Dirac or Majorana. 
The RG evolution of Dirac neutrino masses is studied in 03] . The packages REAP / MPT can 
also be used in this case. 



7.3 Two Higgs Models 

We restrict our discussion to a class of 2HDMs where flavour changing neutral currents 
(FCNCs) are naturally absent [THl [13 [101 • The Yukawa couplings of the theory are given 
by 



i=l 



+ zf rfRFi^^QL0«^ + 4^ uiF»gL0«^} + h.c. , (32) 

where either z^p or z^ has to be zero for each / G {e, u, d, u} in order to ensure the 

absence of FCNCs. In order to generate masses via Yukawa couplings, z^ = 1 for z^ = 
and vice versa. By convention, the right-handed charged leptons always couple to the first 
Higgs, i.e. zi^^ = \,zf^ = 0. 

It is known that in these kind of models there are (at least) two effective neutrino mass 
operators. Furthermore, RG effects are comparatively large, since one has both the tan /? 
enhancement as well as the absence of cancellations due to the SUSY non-renormalization 
theorem. An analytic understanding of the RG effects is more difficult to obtain, since the 
two components of the effective neutrino mass matrix 



= _ ^^(22) (33) 

run differently. Here, more investigations are needed, which are beyond the scope of this 
study. With the REAP package, an extensive numerical analysis is possible. Recently, the 
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RGEs in multi-Higgs models have been derived [HI]. The structure of the /3-functions is 
very similar. 

7.4 Split SUSY 

Let us note that the RGEs for the effective neutrino mass operator in the SM describe the 
running in the framework of 'split supersymmetry' [H21E21 as well, except for a contribution 
to the flavour-trivial part of the RGE (cf. App. lD.3|l . This implies in particular that running 
effects for the mixing angles are suppressed compared to the MSSM (with not too small 
tan/3). The negative g2 contribution to the flavour-trivial part of the RGE gets replaced 
by a positive gl contribution. This effect increases the running of the mass eigenvalues. 

7.5 Other Alternative Sources of Neutrino Masses 

If the dimension 5 neutrino mass operator does not give the leading contribution, possible 
alternative sources to the light neutrino masses can have interesting consequences. Neu- 
trino masses can e.g. emerge from the Kahler potential in supersymmetric theories. It has 
been observed that in this case, large mixing angles can be an infrared fixed point of the 
renormalization group jHHEHj- In the SM, effects of additional dimension 6 operators on 
the running of the dimension 5 neutrino mass operator have been considered in |86j . 

8 Discussion and Conclusions 

We have discussed the running of neutrino masses and leptonic mixing parameters in 
see-saw models involving singlet neutrinos. At energies above the masses of these heavy 
particles, their Yukawa couplings to the left-handed leptons play an important role. As 
they may be of order 1, they can cause significant quantum corrections. We have derived 
approximate renormalization group equations (RGEs) for the mixing angles, CP phases 
and mass eigenvalues. Due to the large number of parameters in the see-saw scenario, the 
details of the running strongly depend on the specific model under consideration. One is 
still able to obtain an extensive analytic understanding of the RG effects. It is instructive 
to compare the RGEs of the physical mixing parameters {ipe} = {9i2,032,023,S,'^i,'^2} 
above the see-saw scales, 

d /^(mfc, phases) j^H.^ v ( i \\ , Mfv v r i \\ f'iA\ 

to those describing the evolution below the see-saw scales. The latter are obtained by 
replacing F^^^iY^, 1^, {il^e}) by Fi{Ye, {ipi}) and F'f^ by zero in Eq. (jHH). Most importantly, 
the structure of the RGEs of the mixing parameters is the same above and below the see- 
saw scales. Hence, there are features common to the evolution above and below. For a 
degenerate spectrum, the first mass quotient in becomes large, yielding strong RG 
effects. There are, however, important differences as well. First, the dimensionless function 
FiiYei {ipi}) vanishes for zero mixing, which is not the case for F'f^ {Yy,Yei{'ipi}). Zero 
mixing angles are hence not stable under the RG in the full see-saw framework. Second, in 
the SM or the MSSM with small tan /5, RG effects are small below the see-saw scales. In 
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contrast, above the entries of Y,y can be of order one and cause important running effects. 
Third, the RGE contains the F^^\Y,y, Y^, {ipe}) term, which describes the radiative rotation 
of Y^Ye in the presence of neutrino Yukawa couphngs Y,y. Finally, between the thresholds, 
there are important effects in non-supersymmetric theories which stem from the different 
scaling of different parts of the effective neutrino mass matrix. 

We listed the leading order RG coefficients for the mixing parameter RGEs in extensive 
tables. Our results allow to obtain a qualitative understanding of generic effects such as the 
influence of the CP phases and that of the absolute neutrino mass scale. For example, non- 
zero phases often damp the running, but some terms in the RGEs are actually enhanced 
by them. A rough quantitative estimate of the size of the RG effects is possible as well. 
Although the change of the mixing angles is quite small for strongly hierarchical masses 
(in the case of a normal hierarchy), it turns out that often it is still comparable to the 
sensitivities of planned oscillation experiments. Therefore, quantum corrections should not 
be neglected in any study of fermion mass models if one aims at theoretical predictions 
whose precision matches that of the experiments. The neutrino mass eigenvalues always 
change significantly due to the RG evolution. This means that a model predicting precisely 
the measured value of Aml^^ = 2.1 ■ 10~^ eV^ at the GUT scale would actually be excluded 
by several standard deviations. Another consequence is a correction to the mass bound 
from thermal leptogenesis. Furthermore, the running of the masses of the singlet neutrinos 
is important for models of resonant leptogenesis. 

In order to obtain precise quantitative results, the complete system of coupled RGEs 
has to be solved. Therefore, one has to resort to numerical calculations. For this purpose, 
we have developed a set of Mathematica packages, which are available at the web page 
http : / / ww w . ph . tu m.de/~r ge/ The package REAP solves the RGEs and thus provides the 
neutrino mass matrix as well as the other parameters such as Yukawa couplings at each 
energy. In models with heavy singlet neutrinos, they are integrated out automatically at 
the corresponding mass thresholds. Thus, the effects of non-degenerate singlet masses, 
which are generally sizable, are correctly taken into account. From the results of REAP", 
MixingPara meterTools, allows to extract the values of the mixing angles, phases and mass 
eigenvalues. 
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Appendix 



A Conventions for Mixing Parameters and Experi- 
mental Data 

A.l Conventions 

Here, we describe our conventions concerning mixing angles and phases. For a general 
unitary matrix we choose the so-called standard-parametrization 

U = diag(e'^% e''^^ ■ \/ ■ diag(e-'^^/^ e-''^^/^ 1) (A.l) 

where 

(C12C13 ■S12C13 Si3e \ 

-C23S12 - S23Sl3Cl2e"^ C23C12 - S23Sl3Si2e"' S23C13 (A. 2) 

S23S12 - C23Sl3Cl2e"' -S23C12 - C23Si3Si2e'^ C23C13 / 

with Cij and Sij defined as cos 9ij and sin6'ij, respectively. 

The MNS mixing matrix Uum is defined to diagonalize the effective neutrino mass 
matrix rriu in the basis where Y^Yg = diag(|/g, ?/^, ?/^), 

Uum^i'UMNS = diag (mi,m2,m3) . (A. 3) 

The mass eigenvalues rrii are positive, and mi < m2 < rris for a normal hierarchy or 
m3 < nil < ^2 for an inverted hierarchy, respectively. For our conventions for extracting 
the mixing parameters from the MNS matrix, we would like to refer the reader to Ref. [20] 
and the documentation of the MixingParameterTools package associated with this study. 

A. 2 Experimental Data 

An overview over the best-fit values and allowed ranges for the neutrino oscillation pa- 
rameters resulting from a global fit to the experimental data [H] is given in Tab. ^| 



Parameter 


Best-fit value 


3(T range 


^12 


33.2° 


28.7° .. 38.1° 


^23 


45.0° 


35.7° .. 55.6° 


^13 


0° 


0° .. 13.1° 




7.9 ■ 10-5 eV^ 


(7.1 .. 8.9) ■ 10-5 eV^ 


l^"^atml 


2.1 ■ 10-3 eV^ 


(1.3 .. 3.2) ■ 10-3 eV^ 



Table 10: Overview of experimental results for neutrino oscillation parameters |41j . 
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B Derivation of the Analytic Formulae 



This appendix contains a couple of technical details relevant for the derivation of the 
analytic formulae discussed in the main part. Our derivation is based on earlier works 
|38| 1^ 1^ . but differs from them by a few steps allowing to express the running of 
the mixing parameters by the mixing parameters themselves rather than mixing matrix 
elements ^2iQj (see also [HZj for real couplings). 

B.l General Strategy 

In an arbitrary basis, one can define unitary matrices f/j, and Ue by 

U^{tf m^{t) Uu{t) = diag(mi(t),m2(t),m3(t)) , (B.4a) 
U,{tyYXit)Ueit) = diag{y',{t),yl{t),yl{t)) , (B.4b) 

where rrij^ is the effective light neutrino mass matrix of Eq. (jH)). The MNS matrix is then 
given by 

Uumit) = Ul{t)U,{t) . (B.5) 
For convenience, we choose to work in a basis, called reference basis in the following, where 

YXito) = diag {ylito),ylito),ylito)) . (B.6) 

Obviously, Ueito) = 1 and UuNsito) = Uu(to). 

Let us now consider the changes caused by changing the renormalization scale according 
to to tQ + At (with At being small). The RGE for rrii, induces a change 

m^{to + At) = m^(to) + ^ [^(^o)^ m^(to) + m^(to) ^(^o) + «(to) m^{to)] 

+ ^((At)2) (B.7) 

with P = {Ce Y^Ye + Cy Y^Yy) in the energy region above the highest see-saw scale. We 
restrict our derivation to this region. As explained in Sec. EJ the results for the MSSM can 

(n) 

also be applied between the see-saw scales after replacing Yy by Yy. However, this is not 
possible in the SM. Due to the change of rriy, 

Uy{to + At) = UM + At Uyito) T + ^((At)2) , (B.8) 

where T is to be calculated below. This relation, however, does not give the full RG change 
of Umns: since also Y^Y^ gets rotated, 

y;ye(to + At) = YXito) + ^ [F\to) YXito) + YX nto) + /(to) YXeito)] 
+ ^((At)2) , (B.9) 

where F = {D^Ye + DyY^u) and / = 2Reae- Hence, C/e(to + At) is different from 
Ueito) = 1 in general, 

Ueito + At) = Ueito) + At Ueito) X + ^((At)^) , (B.IO) 
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with X to be calculated below. 

Using Eq. ()B.5|1 together with Eqs. ()B.8j) and (jB.lOj) . we thus get two contributions to 
the change of the MNS matrix, 

Uumito + At) = UuNsito) + At [UuNsito) T + UuNs{to)] + ^((At)^) . (B.ll) 

We call them the 11^ and the Ue contribution. Following the analysis of 20j, this relation 
allows to derive RGEs for the mixing parameters. 

Before going to the actual calculation, we would like to stress that to derive the mixing 
parameter RGEs, it is useful to work in the reference basis. The resulting equations, 
however, are basis- independent. Of course, if one changes the basis, one needs to transform 
P and F accordingly, which means that the tables in Sec. El and App. O are changed as 
well. 



B.2 RG Corrections Induced by P 

This part of the derivation coincides with the one performed in [20] except for the fact 
that we have to deal with a non-diagonal P. Rather than repeating the analysis of j2()j . 
we just summarize the results: the evolution of U^, is found to be described by 

UlU. = T, (B.12) 

where the entries of T are given by 

levrMmTi,- = -!!^i^^ImP' (B.13a) 

^ rrii + mj ^ ' 

levr^ReT,, = -!I!l±^ReP'. (B.lSb) 
rrii — rrij ■' 

rrii denote the eigenvalues of the effective neutrino mass matrix m,^ (cf. App. IA.1|) . and 
P' = UlPU,. 



B.3 Contribution from the Running of Yg 

Let us now derive the Ue contribution to the RGEs stemming from the fact that Y^Y^ 
changes its structure under the RG. To calculate the corresponding change of the MNS 
matrix, we only need the running of the unitary matrix Ug which diagonalizes Y^Yg. Using 
Eq. (|T^ . it is easy to check that 

IGtt^^YX = F^YX + YXF + '^^^(^eYX- (B.14) 

Plugging this into the inverse of Eq. ()B.4b|) . Y^Y^ = Ue dmg{yl,yf^,yl) f/j =: UeDU^, we 
obtain 

^^{UeDUD = UeDUl + UeDUl + UebUl 

= {F^ U eD Ul + U eD UlF + 2Re aeU eD Ul) . (B.15) 
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Multiplying by from the left and by f/g from the right yields 

UlUeD + D UlUe + D = (F't D + DF' + 2ReaeD), (B.16) 

where F' := FUe- The evolution of Ue can be written as 

^Ue = UeX, (B.17) 

where X is anti-Hermitian. Inserting this relation and using the anti-Hermiticity yields 
i) = (F'^ D + DF' + 2ReaeD) - XD + DX . (B.18) 

By analyzing the off-diagonal parts, we find 

y^X,,-X,,y] = -j^ [iF%,y] + y^Fr] , (B.19) 

where yi = ye etc. For Hermitian F, this can be written as 

y] - yf ^ 

Due to the strong hierarchy of the charged lepton Yukawa couplings, the yi dependent factor 
is approximately ±1. The corresponding equations for the Uu contribution, Eqs. ()B.13|) . 
contain the light neutrino mass eigenvalues, so that a significant enhancement of Tij, the 
analogon of Xij, occurs for quasi-degenerate neutrino masses. In this case, we expect the 
Ue contribution to give only a small correction, unless severe cancellations occur in the U^, 
contribution. However, for a strong normal neutrino mass hierarchy, both contributions 
are generically of the same order of magnitude. The diagonal parts of X, which only 
influence the evolution of the unphysical phases, remain undetermined. 

B.4 Combination of both Contributions 

Inserting Eqs. ()B.12|) and ()B.17|) into Eq. ()B.5|) . we find at t = to in the reference basis 
or 

^MNS ^MNS = T — t/^NS f^MNS ='■ RtX ■ (B.22) 

Note that this is a relation for Uum where both X and T depend on how we split f/jviNS 
into Ue and U,^. Specifically, in an arbitrary basis we have 

uLs Uum = UeT Ul - Ul X U, . (B.23) 

As both sides of the last equation are anti-Hermitian, the derivatives of the mixing param- 
eters are found from the system of linear equations 

+ = Rtx. (B.24) 
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where {^fe} := {9i2, 9is, 623, S, 5^, S^,5r, 1,(^2} ■ The real matrices A^'^^ and S'^'^^ are anti- 
symmetric and symmetric, respectively. Hence, each A^'^) has 3 characteristic elements and 
each S^^-* has 6, so that we can reg ard Eq. ^HMf as a system of 9 linear equations, 

/ Re{RTx)i2 \ 

13 

Re(i?Tx)23 

Im(i?Tx)ii 

Im(i?Tx)i2 . (B.25) 
Im(i?Tx)i3 
Im(i?Tx)22 

Im(i?Tx)23 
y lm{RTx)33 j 
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12 


. 4(8) 

12 


\ 


4(1) . 
lo 


lo 




4(1) . 

23 


23 




s['^ ■ 

11 


11 




0(1) 

*_? 1 


0(9) 
12 




q(1) 
^13 


0(9) 

^13 




^22 


0(9) 
^22 




'-'23 


0(9) 
^23 






0(9) 
■ "^33 


/ 




= : V 



V can be split into two parts, 

V = Vt + Vx , 



(B.26) 



where vt is built from T and vx is built from —Ul^^^X Uum- In particular, each ^k, for 
instance 612, is the sum of two contributions, one from T (i.e. from the running of m,y) and 
one from X (i.e. from the running of Ye). 

B.5 Comment: 'Unphysical' Phases 

The RGEs in the full theory contain the entries of P. However, the phases appearing in 
the off-diagonal elements of P are not basis-independent, rather they can be changed by a 
transformation using the 'unphysical' phases 6e, 5^ and 5r only. To see this, let us perform 
(in the basis where Y^Yf, is diagonal) a transformation /T, 



K 



where K = disig{e''''\ e''^\ e'^^ 



Ken, (B.27) 
) is a diagonal phase matrix. Y^Y^, is invariant under this 



transformation, yet it changes the effective neutrino mass matrix according to 
K* rrijj . 

Hence, also f/jviNS gets changed under this transformation, 

K 

UmNS — ^ -^f^MNS 5 



(B.28) 



(B.29) 



i.e. K affects the phases 6e, 5^ and 5-r in the standard parametrization (jA.lj) . Furthermore, 
it rotates the phases of the off-diagonal entries of Y^Y^j as 



K 



KYjY.K^ . 



(B.30) 
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This shows that one has to specify both the phases 5e, S^, Sr and the arguments of the 
off-diagonal entries of YjY,y, as one set of parameters can be traded for the other. In other 
words, two theories with equal P but different phases 5/ are not equivalent. In the main 
text, we use the convention 

Se = 6^ = 5r = . (B.31) 

As a technical comment, we would like to mention that in order to diagonalize a gen- 
eral neutrino mass matrix m^, the parameters 6^, 6^ and 6r are needed. Only after the 
transformation with K = diag(e~^'^% e"^'''', e"^''^), one can write the MNS matrix without 
Se, Sfj, and Sr- The step of going to the basis where 6e = Sfj, = 6r = has often not been 
mentioned explicitly in the literature. In this context, we would like to comment that, of 
course, 6e, Sfj, and 6r are subject to quantum corrections with their RGEs depending on the 
physical parameters. 6e has a term proportional to l/Am^^i whereas 5^ and 6r are both 
proportional to 1/Am1^^}^ 



C RGE Coefficients 

In the following, we show the RGEs for the lepton mixing parameters obtained from 
the derivation discussed above. We give the first order of the expansion in the small 
CHOOZ angle 6'i3. We furthermore use the abbreviation ( for the ratio of the mass 
squared differences, cf. Eq. (fTHjl . 

The results are presented in the form of tables which list the coefficients of Pjg = 
{CeYjYe + CyY^Yy) fg in the RGEs. Thus, if only a single element of P is dominant, the 
derivatives of the mixing parameters are found from the corresponding rows in the tables. 
Of course, if several entries of Pjg are relevant, their contributions simply add up. While 
the complete RGEs are basis-independent, the table entries do depend on the choice of the 
basis, since P is basis-dependent. We use the basis where Y^. is diagonal and where the 
unphysical phases in the MNS matrix are zero. 



^®The corresponding formulae below the see-saw scales can be obtained from the web page 
"http : //www . physik . tu-muenchen . de/~mratz/Analyt icFormulae/ There, the RG evolution of the 5f 
phases depends on the physical parameters, but the RGEs of the physical parameters are independent of 
the 5f phases. 
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Q 
Q 



_j_ 


1 < / tg 1 lL\t!. -'- \ 


13 — 


Am^, (1+C) 




|m3±m2e''^2|2 


23 — 





|m2e''^2 ilimie''^i p 

sol 



<Sl3 
<S23 
<Sl2 



mi?Tt3 sin 
m2m3 sin 1^2 
mim.2 sin(</Pi-i/;2) 



Am" 



ol 



._|_ (m^+mg) cos 5it2mim3 cos((5—ipi) 

•^13 = Am^,^(l+C) 

(m|+m|) cos (5it2m2»n3 cos((5— 1^2) 



•^23 



Ami 



E 



_|_ (mj+mg) sin (5ib2mim3 sin(<5—(pi ) 

13 = Am^,^(l+C) 

(^m^+mg^ sin <5it2m2m3 sin((5— (/P2) 

23 ^ a;^?;;:^ 



/^12 
'-'13 



/^23 
*-'13 



/^12 
'-'23 



ri3 

'-'23 



Am"(i+c) + ^) sin ((^1 - V22) - Crn3 sin (25 - (^i)] 



sol 



[nil sin {26 - + (1 + C) "^■2 sin (p2] 



[mi sin (v^i - (p2) - (1123 sin (25 - V52)] 



sol 



Am2,^(l+C) 



[mi sin (/?i + (1 + C) m2 sin (25 - (P2)] 



V2 



'atmV 



[mi COS (5 — v^i) — (1 + C) m2 cos (5 — (p2)] sin 5 



Am2,„,(l+C) 



[mi cos (25 - ipi) - (1 + C) m2 cos (25 - (P2) + ("^3] 



Table 11: Definition of tlie abbreviations used in Tabs . 1^ and 1 1 21 [TH 





327r2 012 


647r2 ^13 


32712^23 


Pll 

P22 

P33 

ReP2i 
RePai 
ReP32 
ImP2i 
ImPai 
ImP32 


Q^2 sin 2^12 
sin 26*12023 
-Q^2 sin 26*12^23 

2Q+2C0S26^12C23 
-2Q^2 cos 26*12^23 

Qf2 sin 26*12 sin 26*23 

45i2C23 
— 45i2S23 






(^^3 - ^^3) sin 26*12 sin 26*23 
-(^^3-^13) sin 2^12 sin 2^23 
4(^+022 + ^+522)523 
4(^+022 + ^+522) C23 
2 (^+3 - ^+3) sin 2^12 cos 2^23 

4 (i3i-3c22 + S2-3S22) S23 
4 (^r3C?2 + I323SI2) C23 

2 (-62^3 - i3i"3) sin 2^12 




(Q^3C?2+ Q^3S?2) sin 2^23 
-(Q2V?2 + 2l35?2) sin 2^23 
(2^3- 2^3) sin 26^12523 

(2ii- 2^3) sin2^i2C23 

2(Q2V?2 + 2^3^?2) COS 2^23 

2 (^23 - ^13) sin 26*12^23 
2(^23 -5i3)sin26*i2C23 

4(523C22+513S22) 



Table 12: Coefficients of Pfg in the RGEs of the mixing angles 0ij in the limit ^13 0. The abbreviations 
Afj, Bfp Sij and Q+ depend on the mass eigenvalues and phases only, and enhance the running for a 
degenerate mass spectrum, since they are of the form /i^ (m^ , TOj , phases) /(m^ — m|). They are listed in 
Tab.ini 
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647r25(-i) 


Pll 





P22 


- (i3+3 - B+) sin 2^12 sin 2^23 


P33 


(-6^3 - Bt^) sin 2^12 sin 2^23 


RePai 


-4 {Bt,cl2 + S2+3^?2) ^23 


RePsi 


-4 (i3+3C?2 + S2V?2) C23 


ReP32 


-2 (-6^3 - Bts) sin 2^12 cos 2^23 


ImP2i 


4 (^i3Ci2 + -423-512) -523 


ImPgi 


4 (^r3C?2 + ^235?2) C23 


ImPgs 


2 (^23-^4:3) sin 2^12 





647r25(o) 


Pll 
-P22 
-P33 

RePai 
RePai 
ReP32 
ImP2i 
ImPsi 
ImP32 


-8 iiCa + 5i2 - ^23) C?2 + (ClI + 5i2 - 5i3) S?2) 
8 (((5i2 - 523) 4, + Ci.^3) C?2 + ((5i2 - 5i3) 0^3 + C^sl,) sl,) 
8 ((^13^23 + ('^12 - '^23) S23) C^2 + (C2IC23 + {S12 - ^13) S23) S12) 

-165i2C23 cot 2^12 + 4 (21)1023 + (^23 - ^13) S23 tan 6^23) sin 2^i2 
165i2S23 cot 2^12 - 4 {2V1S2Z + (^23 - Si^) C23 cot ^23) sin 26^12 
-16 (523C?2 + Si^isl^) COS 2^23 cot 2^23 - 8 {elicit + Cllsl^) sin 2^23 

-8Qr2C23 CSC 2^12 - 2 (2P2C23 + {Q23 - Qvi) COS 2^23 SCC ^23) sin 2^12 

8Qr2S23 CSC 2^12 + 2 (2P2S23 - (Q23 - Qu) COS 2^23 CSC ^23) sin 2^12 

-8 (Q2"3C?2 + Qr3S?2) cot 2^23 



Table 13: Coefficients of Pfg in the derivative of the Dirac CP phase. The complete RGE is given by 
5 = 6'f3^(5(~i) +(5^°) + ^(6*13). The abbreviations Af^, B^j, Qf^, C^j and depend on the mass eigenvalues 
and phases only, and are listed in Tab. 1111 
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167r^0i 


Pll 


-45i2c22 


P22 


45i2C?2C23 - 4 (523C12 + Si^sl^) COS 26*23 


-fss 


45i2C?2S23 + 4 (523C?2 + 5i3S?2) COS 2^23 


ReP2i 


-45i2C23 COS 2^12 cot 9i2 - 2 (523 " Si'i) COS 2^23 sec ^23 sm 2^12 


RePgi 


45i2S23 COS 2^12 cot 9i2 - 2 (^23 " Si'i) cos 26*23 CSC (923 sin 2^12 


ReP32 


-8 (523C12 + S12.SI2) COS 26*23 cot 26*23 - 45i2Ci2 sin 26*23 


ImP2i 


-2Q5"2C23 cot 6*12 - {Q22, - Q12) COS 26*23 sec 6*23 sin 26*12 


ImPgi 


2Qr2«23 cot 9i2 - (Q23 - QT^) COS 2^23 CSC ^23 sin 2^12 


ImP32 


-4 (Q23CL + 2r3^?2) cot 2^23 






167r^<^2 


Pll 


-45i2S?2 


-P22 


ASi2cl^sl^ - 4 (>S23C?2 + ^5135^2) COS 2^23 


-P33 


45i2si3S?2 + 4 («523C?2 + '5l3S?2) COS 26*23 


ReP2i 


-45i2C23 COS 2^12 tan 9i2 - 2 (^23 - ^13) cos 2^23 sec ^23 sin 2^12 


RePsi 


45i2S23 COS 2^12 tan 6*12 - 2 (523 - Si^) cos 26*23 esc 6*23 sin 26*i2 


ReP32 


-8 (523C12 + 513S12) cos 26*23 cot 26*23 - 4^125^2 sin 26*23 


ImP2i 


-2212^23 tan 6*12 - (Q23 - Qiz) cos 26*23 sec 6*23 sin 26*12 


ImPai 


2Qr2*23 tan 6*12 - (Q23 " Qr3) cos 26*23 esc 6*23 sin 26*12 


ImPgs 


-4 (Q23C?2 + Qr3S?2) cot 2^23 



Table 14: Coefficients of Pfg in the RGEs of the Majorana phases for ^13 = 0. 
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D RGEs for See-Saw Models 



In order to calculate the RG evolution of the effective neutrino mass matrix, the RGEs for 
all the parameters of the theory have to be solved simultaneously. We therefore summarize 
the RGEs for the minimal see-saw extensions of the SM, of the class of 2HDMs described in 
Sec. O and of the MSSM. We list the MS 1-loop results in the SM and 2HDM, as well as 
the 2-loop RGEs for the effective neutrino mass operator, the singlet mass matrix and the 
Yukawa couplings in the MSSM. For further RGEs and references, see e.g. f88. 89. 90, 9J]- 
We use the notation defined in Sec. [21 In particular, a superscript (n) denotes a quantity 
between the nth and the (ra + l)th mass threshold. The RGEs for the SM, 2HDM or MSSM 
without singlet neutrinos can be recovered by setting the neutrino Yukawa coupling to zero. 
In the full theories above the highest see-saw scale, the superscript (n) has to be omitted. 

The RGEs for the gauge couplings are well-known and not affected by the additional 
singlets at 1-loop order. They are given by 



Ifi 2 

lovr jj,— — 
d/i 



bAg\ 



(D.32) 



with f6sU(3)c) ^SU(2)l) ^U(1)y' 
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-7,-f ,fi) in the SM, (-7, -3, f ) in the 2HDMs and 



3, 1, ^) in the MSSM. For U(1)y, we use GUT charge normalization. 



D.l The RGEs in the Extended SM 

In the SM extended by singlet neutrinos, the RG evolution is governed by the /5-functions 

Emm 



(n) 
(n) 



An) 



(n) 



+ 2 Ti{Y^Ye)'K + 2 Tr (PfY^YK + 6 TiiY^Y^fK 
+ 6TriYjY,f^-3gf^+)il^, 

An) (n) (n) (n) ,{n) (n) 

to) M + M (ni;t) , 

(n) ( 3 An) (n) 3 An) (n) 

Y.\^^{yX) - ^(YX) + Tr (yX) + Tr(F;ye) 
+ 3 Tt{YXu) + 3 Tr(F;F,) - ^gj - ^-glj , 

-'^p ^ 2^ 2 4 4 



+ Tr 



YXe + YjYu + 3 yJy, + 3 YXu 



Y.[IyX-IyXu-\91-191-891 



+ Tr 



(n)(n) 



YXe + YXu + 3 YjYd + 3 YXu 



(D.33a) 
(D.33b) 

(D.33c) 



(D.33d) 



;D.33e) 
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Y /-yty _ ^yty II 2 _ 9 2 _ o 2 
-'u ^ 2 " " 2 20 4 

(n) (n) 

F^Fe + + 3 FjFd + 3 yjK 



+ Tr 



(D.33f) 



+ 4A Tr 
- 8 Tr 



(n) (n) 

Y}Y, + FjF, + 3 Y}Yd + 3 



(D.33g) 



F^Fe nt^e + YX YX + 3 Y^Yd YjVd + 3 



We use the convention that the Higgs self-interaction term in the Lagrangian is — 



D.2 The RGEs in Extended 2HDMs 

Here, we hst the /3-functions for the class of 2HDMs described in Sec. 17.31 j^ lH^ ISTj. The 
coefficients z^^^ determine which fermion couples to which Higgs, cf. Eq. 



(ii) 



(") 

167r2/3M 

J") 



2 Tr(Fjn) 



(1 _ 25a) W''\yX) + {YXfn^^^^ + 
+ 6a 2 Tt{YX) + 6 Tt{Y^Y^) + zf 6 Tr(r;F, 
+ AiK^") + ^aA^K^^^) + ^isAsK^^^) - 3^2k(**) , 

(n) (n) (n) (n) (n) -Ji 

(nF;)M + M(F.i;t) , 



+ Y/^ Tr 



^aV-^i; + MIH + ^zf yJy, + 3z« rjy; 



IQ-'Py. = Y,[^YX+{l-2z^;^)YX-y,-y2 



+ Tr 



YX + il^E + 34^) YjYd + 34'^ YXu 



+ E4^^Tr 



i=l 



An) In) 



dXYe + z^yXu + ^YlYd + 34)i;^y; 



(D.34a) 
(D.34b) 



(D.34c) 
(D.34d) 



(D.34e) 
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(n) 



167r^ 



+ Tr 



1=1 



(D.34f) 



For the parameters of the Higgs interaction Lagrangian, the /5-functions are [92j (Note that 
we use different conventions for the renormahzable Higgs couphngs, as specified in ^7].) 



(n) 



167r^ (3, 



Ai 



(n) 



167r^/3A3 



(n) 



6A? + 8A2 + 6A3A4 + A^ - 3Ai [?>gl + '^gl^ + 3^2^ + ^ {^-gl + 
+ 4Ai Tr {y}Y, + z(^) y;t?>;^ + 341) yjy, + 3^1) yjy^^ 

((n) (n) (n) (i) /-, % , , \ 

(D.35a) 

Q\l + 8A2 + 6A3A4 + A^ - 3A2 (^gl + + ^gt + ^ + gl 

(n)ln) 



+ 4A2 Tr (^42) rj-yt + 3^'^ yJy^ + 34') Fjr^ 

- 8 Tr y;t5>^5>yy;^ + 342) yJy,Y}Y, + 3z^^^ Y^Y^Y^Y^^ , (D.35b) 
(Ai + A2) (3A3 + A4) + 4X1 + 2A^ + - 3A3 (^3gl + ^gl^ + ^g^ 

+ Y^gt - ^glgl + 4A3 Tr [yX + yX + 3i;V, + ^Yuj 

- 4 Tr YX Yit + 3 (z^^^zi^) + zf z^^^) Y^ Y^u^ , (D.35c) 



IQn' P^, = 2 (Ai + A2) A4 + 4 (2A3 + A4) A4 + 8A^ - 3A4 3^72' + -gf 



(n) 



(n) (i) 

+ 4A4 Tr YX + YXu + XYd + 3KjK 



' 2 2 

r^i^2 



+ 4 Tr ( 4^)Fiye yVy;^ + 3 ( + 42)41)) yjy^ yXu 

'3 



(D.35d) 



A, 



Ai + A2 + 8A3 + I2A4 - 6 ( -gl + ?>gl 



(n) (n) 

2 Tr F^Fe + YXv + 31^;^, + 3rjF„ 



(D.35e) 



D.3 Split Supersymmetry 

The /^-functions for the renormahzable couphngs in the framework of spht SUSY are 
hsted in Ref. (HB]- The diagrams contributing to the RGE of the effective neutrino mass 



41 



Or -rO- 



B 



Figure 9: Additional diagrams contributing to the wavefunction renormalization of the Higgs in spht 
SUSY. The Higgsino is denoted by </>, and B and ]¥"■ represent the Bino and the Winos. 

operator are those relevant in the SM, amended by two diagrams involving Higgsinos and 
gauginos (cf. Fig. IH)). These diagrams contribute to the flavour-trivial part of the RGE. At 
1 loop, we obtain for the divergent parts of the renormalization constants in dimensional 
regularization and in the MS scheme 



Sk 1 



1 



167r2 



2 Tr (SY^ + 3YjY, + YX) + ^(^i - 1)^7? + ^(6 - 1^ 



167r2 



2k {YX) + 2{YXf>^ 



K 



where and ^2 are the gauge parameters in gauge. Zi := 1 + 6Zi^i 7 + • • • (e := 4 — rf) 
are wavefunction renormalization constants, and 6n := 6k^i ^ + . . . is defined via the 
counterterm for the dimension 5 operator. 



■ 6Kf, . 0) (il . 0) + h. 



Using the method described in we then find the 1-loop /5-function 



+ 



X + -gf + 2 Tr (syX + + YX 



(D.37) 



Clearly, the term involving the gauge couplings in the flavour-diagonal part differs from 
the SM case. 



D.4 The RGEs in the MSSM Extended by Heavy Singlets 

f(n) (") (") 1 

We give the 2-loop RGEs for the quantities Q E < k,M, Y^, Y^, Y^, Yg } in the form 



(D.3^ 
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The 1-loop parts are given by HZ] 



(4>r)^/3i'> 



(Ftn)-'s+'S(Fty.)+(?;t'fy'K>+%>(yyF:) 

+ 2Tr (y;;i'')'A'+ 6Tr(y;ty„)S'- ^sfS- 6g| 

5 

/(") ("), N (n) (") (i), \ T 

2(ni;t)M + 2M(r,r;) , 



(n) 

K , 




(D.dyaj 
(D.39b) 


3 o 




(D.39c) 


7 

' 2 


- ^gl - 


^P. \ 
2 1 

(D.39d) 




- ^gl - 




9 2 

5^^- 




(D.39e) 


3^72'} • 


(D.39f) 



The results for the 2-loop parts, which are an extension of the usual 2-loop /3-functions 
for the MSSM jHS!, are jM! 



(n) 



6Tr(F„r;r,Fj) - 18Tr(r„Fj"y„yJ") - 2Ti{YXyeY}) 



(n) (n) (n) (n) S 007 1 Q 

6 Tr(FjnFjn) + -^7? Tr(F„tr„) + ^2gl Ti{Y^Y^) + — + -glgl + 15^?: 
5 25 5 



(n) 



(n) (n) (n) (n) 



2(F;F,y;F,)^ + 2(yjnyjy,)^ + Tr(ni;t) + 3Tr(F„rJ) {yX) 



1) ("), 



+ ( -^1 + Tr(i;i;t) + 3Tr(r,Fi) ) {YJY^f 



(n) 



(») 
- K 



(n) (n) (n) (n) / (n) (n) \ (n) (n) 

2r;rei;^re + 2rjnrjr, + Tr(r,r;) + 3 Tr(r„y;t) 



6 



+ ( --gl + Tr(y;r;) + 3 Tr(y,rJ) ) y;ty; 



(47r)^/3iJ) = M 



- 2r;FjF;Fj - 2y:yJy:yJ - 6y:yJ Tt{YX) 



+ 



(n) (n) (n)(n) 5 {n)(n) 

- 2Y:yJ Tt{YX) + -gl Y:yJ + 6gl YJYJ 

(n) (n) (n) (n) (n) (n) (n) (n) (n) (n) (n) (n) 

- 2Y,yXyJ - 2ny;tnF; _ Tt{yX) - 2nrJ Tr(r,y;t) 



6 „(")("), „(")("), 

+ -glYX + QglYX 



(n) 

M 



(D.40a) 



(D.40b) 
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(n), (n) I (n) (n) (n) (n) (n) (n) 

(47r)^/3^f = yul- 2YXy}ye - 'lYlYeYlY. - ^Y^Y^Y, - ?>YX Tr(ldrJ) 

(n) (n) (n) (n) (n) (n) (n) (n) 

- y^Fe Tr(Fei;^) - ?,yX Tr(r.r;) - grjy; ti{YuYI) - Tr(y;i;ty^yj) 

(n) (n) (n) (n) 6 6 

- 3 Tr(yJnFjy;) - 3 Tr(y;rir,i;t) - 9 Tr(y;y;ty;y;t) + y}y, + -^7? rjy. 

+ ^glY^l + ^9! TV(F„tyj + 16^3^ Tr(i;tF„) + ^gf + ^glgl + ^g\ } , 

(D.40c) 

(47r)^^J.f = Ya^- ^yIy.yIy, - 2yXyJy, - 2yXyX " 9 Tr(y,rir,r]) 

(n) (n) 

- 3Tr(F,Fjr„Fi) - 3 Tr(rer;F,i;t) - TriY^Y^Y^^) - 9Y^Y, Tt{Y,yJ) 

(n) (n) 

- 3YjY, Tr(F,y;t) - Y^ MY^) - 3YX MYuY^) + 6gl yJy, 

+ \glY}Y, + ^glY^Y^ - ^gf TriY^Y,) + ^gf Tt{YX) + ^QqI MyJy,) 
+ ^-^9t + 9l9l + §9t + l9hl + Sgh! - fgt } , (D.40d) 
= Y^ [-2Y^Y,YjY, - 2YJy,YX - AY^Y^ - SYJy,Tt{Y,Y^) 

(n) (n) 

- YjY, Tr(i;F,t) _ QY^Y^ Tr(F„i;t) _ gyjy^ Tr(ni;t) 

(n) (n) (n) (n) (n) (n) 

- 3Tr(r„r]r,rJ) -9Tt{Y^Y^Y^Y^) - Tt{y,yXyJ) -3Tt{Y,YJY,YJ) 

+ ^gf YlY, + yi YlY^ + ^gl Y^ + Tr(i;ty;) + iQgl Tv{Y^Y.,) 

2743 4 9 9 15 4 136 o o „ o o 16 , 1 /-r^ x 

+ + ^7?^?2 + Y^?2 + ^^?i^73 + ^ahl - -^gt | , (D.40e) 

(47r)^/5^;) = Ye\- AYXY^Y, - 2YXyX ' '^YJY^YX - ^Y^ Tt{Y,yJ) 

(n) (n) (n) (n) (n) (n) 

- 3YX MYeYj) - YX Tr(F.F;) - 3YX Tr(r„Fi) - 9 Tr(F,r;F,Fj) 

- 3 iTiY.Y^Y^Yi) - 3 TiiXXYeY}) - TiiY^YX) + ^qI ^Y^Y^) 

+ Qgl YX - ^riYjY,) + IQgl Tr{YjY,) + ^^gf + ^g^gj + ^g^ } . 

(D.40f) 

The 2-loop /5-functions for the gauge couphngs in the presence of Y"^ can be found in |29j . 
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